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involves the propagation of massive quasiparticles. The appearance of complex poles is interpreted 
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The Nambu- Jona-Lasinio model (NJL) has been vastly 
considered for studying nonperturbative aspects of QCD. 
Nowadays, it is mainly used to explore finite tempera- 
ture and density effects in the frame of the mean-field 
approximation [lJ-Q . One of the big challenges in QCD 
is to understand the confinement mechanism and the dy- 
namics behind confinement. Perturbative QCD cannot 
describe confinement and, although lattice QCD is able 
to reproduce successfully hadron properties, like masses 
and coupling constants Q, it has problems when deal- 
ing with finite baryon chemical potential (the sign prob- 
lem). However, there are effective models which include 
explicitly confinement, as, for example, different versions 
of the bag model [fljllj Dyson- Schwinger models [9rJl3j|. 
or the Polyakov loop effective action coupled to Dyson- 
Schwinger or NJL models [I3 - [r7l ]. 

The nonlocal NJL model (nNJL) is another attempt in 
this direction fTsM2l1 | . When the gluon degrees of freedom 
are integrated out in the QCD action, a nonlocal quark 
action emerges and confinement should be hidden there. 
The idea of the nNJL approach is to incorporate nonlocal 
vertices through the presence of appropriate regulators. 

Since the NJL model is nonrenormalizable, a momen- 
tum cutoff is needed in order to handle the UV diver- 
gences. The applicability of the model is, therefore, 
restricted to energy scales below the cutoff. Nonlocal 
extensions of the NJL model are designed to regular- 
ize the model in such a way that UV divergences are 
controlled, internal symmetries are preserved, and quark 
confinement is incorporated. The nNJL model has been 
extended to a finite temperature and density scenario 
[23 - l32j through the usual Matsubara formalism [H, [34| . 
Here, nevertheless, the exact summation of Matsubara 
frequencies turns out to be cumbersome, due to the com- 
plicated shape of the regulators. In most cases, it is nec- 
essary to cut the series at some order. 

The idea of this work is to develop the finite temper- 
ature real time formalism for the nNJL model. In this 
way, we are able to calculate temperature corrections, 



providing a physical picture in terms of quasiparticles. 
On the one hand, those states with real masses can prop- 
agate freely in the deconfined phase. On the other hand, 
the existence of complex poles of the propagator in the 
confined phase produces a strong damping avoiding the 
propagation of such states. 

Real time formalisms, as thermo-field dynamics or 
the Schwinger-Keldysh formalism, can be constructed 
through the analytic continuation of the Euclidean action 
[33M39| . Those formalisms double the number of degrees 
of freedom, providing the appearance of a 2 x 2 matrix 
propagator SV,-. In any version of the NJL model in the 
mean-field approximation, the gap equation corresponds 
to a one-loop self-consistent relation, so we only need to 
find the S\i component 0, The other matrix prop- 
agator components start to appear at the two-loop level. 
In the construction of the Sn propagator, the main ingre- 
dient is given by the spectral density function. However, 
the construction of the spectral density function is not a 
simple task in a nonlocal frame, especially when dealing 
with a nontrivial analytic structure of the propagator. 

In this article, we will develop the real time for- 
malism for the Gaussian regulator in the nNJL model, 
which can easily be extended to other kind of regulators. 
As a result, we will get thermal propagators describing 
quasiparticles with temperature- and chemical-potential- 
dependent masses. 

This paper is organized as follows: In Sec. U we intro- 
duce the nNJL model in the mean-field approximation at 
finite temperature and chemical potential. Section |TT] is 
devoted to the discussion of the analytical structure of 
the nNJL model in the mean-field approximation with a 
Gaussian regulator. In Sec. IIII1 the real time formalism 
is presented for the general case, constructing, then, in 
Sec. IIV1 the corresponding real time propagator Su for 
the case of a Gaussian regulator. We explain how to ob- 
tain the gap equation and the thermodynamical potential 
in Sec. IVland Sec. I VII respectively. Finally, in Sec. I VIII 
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we summarize our conclusions. 



I. nNJL MODEL IN THE MEAN-FIELD 
APPROXIMATION AT FINITE TEMPERATURE 
AND CHEMICAL POTENTIAL 

Nonlocal models with separable interactions have been 
considered several times. In the context of NJL, confine- 
ment is introduced through the inclusion of interacting 
nonlocal currents, which are extended regularized ver- 
sions of the usual local currents. 

The SU(2) f ® SU(3) C nNJL Euclidean Lagrangian is 
given by 

G 

Ce = ip(-i0 + m)ip - -^ja(x)j a (x), (1) 

m being the current quark mass and j a (x) the nonlocal 
quark currents. Here, we are using the metric — 
diag(l, 1,1,1). 

There are two schemes on the market to introduce non- 
local effects in terms of extended currents [2l| . We will 
use the one based on instanton liquid models: 



ja(x) 



d 4 y d 4 z r(y - x)r(z - x)ip(y)T a '4>(z), (2) 



where r(x — y) is the nonlocal regulator and T a = 
(1, 1757"), r being the Pauli matrices. 

Since we want to deal with mesonic degrees of freedom, 
we will follow the standard bosonization procedure. This 
is realized through the introduction of auxiliary scalar 
and pseudoscalar fields. By integrating out the quark 
fields, an equivalent partition function, in terms of only 
bosonic degrees of freedom, is obtained: 



Z = I DaD 3 TT e- T[a ^ a 



(3) 



We proceed within the mean-field approximation, keep- 
ing the mean values of the boson fields and neglecting 
their fluctuations. In this way, the partition function in 
Eq. ([3]) turns out to be Z « Zmf = e~ rMF , where the 
mean-field effective action in Euclidean momentum space 
is given by 



MF 



Va 



a 

2G 



2N r 



(2^ trln5 * 



Gte) 



(4) 



Qe = {q, 94) being the four-momentum in Euclidean 
space and V4 the four-dimensional volume. The trace 
acts on the Euclidean Dirac matrices 75 = (7,170) of 
the effective quark propagator 



S E (qi 



-^ + £( gg 2 ) 
q E 2 + X 2 (q E 2 )' 



(5) 



where the running mass includes the Lorentz invariant 
regulator contribution in Euclidean momentum-space 



S (<?-E 2 ) =m + ar 2 (q E 2 ), 



(6) 



a being the mean-field value of the scalar field. The pseu- 
doscalar mean-field value is zero due to isospin symmetry. 
Chiral symmetry, on the other hand, is explicitly broken 
through the current quark masses and spontaneously bro- 
ken by a nonvanishing chiral condensate value. 

Finite temperature (T) and chemical potential (//) ef- 
fects are introduced in the standard way through the 
Matsubara formalism [33j]. As a result, in Eq. (j4j, the 
four-dimensional volume, the momentum, and the inte- 
gral in momentum turn out to be 



Va -> V/T, 

<74 ->■ -q n , 

dqj 
2ir 



(7) 

(8) 

(9) 



where q n includes the Matsubara frequencies and the 
chemical potential 



q n = (2n + l)nT + ifx. 



(10) 



The grand canonical thermodynamical potential is 
given by n MF {a,T 7f i) = {T/V)T MF {a,T^) g|. The 
value of a is obtained then through the solutions of 
9Qmf/9ct = 0, giving rise to the gap equation 



j0r 2 (,, 2 ,+<, 2 )trS £ ta,-,„). (11) 



This equation, which involves directly the propagator, 
will be our starting point in order to apply the real time 
formalism. 



II. POLES OF THE EFFECTIVE QUARK 
PROPAGATOR WITH A GAUSSIAN 
REGULATOR IN MINKOWSKI SPACE 

Among the different kind of existent covariant regula- 
tors, we will use here the Gaussian one due to its simplic- 
ity [18| |. In Euclidean momentum space, the regulator is 
given by 



r(q E 2 ) 



exp 



QE 
2A 2 



(12) 



where A is a free parameter that has to be chosen from 
phenomenological considerations. Thus, the effective 
model depends on three parameters: the current quark 
mass m, the effective coupling G, and the scale parameter 
A, the last one being associated to the cutoff in the usual 
NJL model by setting the regulator as r = 8(A 2 — q 2 ). 
These parameters are fixed in order to get the physical 
values of the pion mass m n = 139 MeV, the pion de- 
cay constant f„ — 92.4 MeV, and the chiral condensate 
-(qq} 1/3 ~ 200-260 MeV Following [H, we will use 
two sets of parameters. Set I is given by (qq) = —(200 
MeV) 3 , m = 10.5 MeV, G = 50 GeV" 1 and A = 627 
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FIG. 1. Schematic description of the poles of the propagator 
in Minkowski space. Mi and M2 correspond to the deconfined 
masses, M c is the critical mass and Mq and To are the mass 
and the decay width in the confined phase, respectively. 



MeV. Set II is given by (qq) = -(220 MeV) 3 



7.7 



MeV, G = 30 and A = 760 MeV. 

To start our analysis in the real time formalism, let us 
explore more in detail the analytic structure of the prop- 
agator by setting q^ = iqo. The Euclidean propagator in 
Eq. ([5]) turns out to be Se — > iSo, where the propagator 
in Minkowski space is defined as 



So(q) 



: t+n-1 2 ) 

'q 2 -E 2 (-q 2 y 



(13) 



and where now q = (qa,q) and 7 = (70,7), with the 
metric g^ = diag(l, — 1, — 1, — 1). This propagator, as 
shown in Fig. [1] presents three different kind of poles 
which depend on a. For Low values of a, there are two 
real poles associated to real masses that we denote as 
Mi and M%. In this case we will say that the system is 
in the deconfined phase since we have freely propagating 
quasiparticle states. As a grows up to a certain critical 
value, the two masses join into a single real mass M c . For 
higher values of a, this critical mass splits into complex 
poles. There exist an infinite number of those poles for 
the Gaussian regulator. Nevertheless, it can be shown 
that the only relevant poles correspond to the first pair 
[HI], where the real part is associated with a mass Mo 
and the imaginary part is related to a decay width factor 
Tq. The other complex poles involved have considerable 
higher values for the masses and decay widths, and these 
masses are not continuously connected with the critical 
mass M c . The appearance of c omp lex poles is interpreted 
as a signal of confinement [III Ell! since the correspond- 
ing quasiparticles do not propagate freely. As soon as the 
poles becomes complex, we will be in the confined phase. 

The poles of the propagator in Eq. (fTB")) are plotted 
in Fig. [2] for both sets of parameters as a function of 
the self-energy at zero 4th-momentum E(0) = m + a. 
When E(0) = S c , there is only one positive real pole 
q 2 = M 2 , which is determined from the condition d q 2 [q 2 — 
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50 100 150 200 250 300 350 400 
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FIG. 2. Poles of the propagator as a function of the self- 
energy at zero 4th-momentum for set I (upper graph) and set 
II (lower graph). Defining the complex poles as q 2 = Mq ± 
iMoTo, the solid, dashed, and dot-dashed lines correspond 
to the masses, and the dotted line corresponds to the decay 
width. 



E 2 (-g 2 )] =0, giving 



E c = m + (M c - m)e- M = /A 
-(m + \fn 



Mr = 



2A 2 



(14) 
(15) 



The corresponding values are E c — 273 MeV for set I 
and E c = 329 MeV for set II. The previous analysis can 
be easily verified by plotting q 2 — £ 2 (— q 2 ) as a function 
ofq 2 . 

The values of S(0) at zero temperature and chemical 
potential are E = 350 MeV for set I and £ = 300 MeV 
for set II, respectively. The tendency for the a parame- 
ter is to decrease as temperature and chemical potential 
increase. Then, if at zero temperature and chemical po- 
tential the corresponding poles are real, they will still be 
real at finite T and \x. On the other hand, if the poles 
are complex at zero temperature and chemical potential, 
they will become real at some finite values of T and [i. 
As £ > E c in set I and £ < E c in set II, they are called 
confining and nonconfining sets, respectively. 

Now we proceed to extend this formalism to a finite 
temperature and density scenario in the frame of real 
time formalism. 
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III. REAL TIME FORMALISM 

It is "well-known that when formulating quantum held 
theory at hnite temperature in the real time formalism, 
we have to double the number of degrees of freedom (33l — 
l39j . The new fields that appear are called thermal ghosts. 
As a consequence, the new thermal propagators are given 
by 2 x 2 matrices Sy. However, for one- loop calculations, 
we only need the term Sn. The other components of 
the thermal matrix propagator participate only in higher 
loop calculations, since the thermal ghosts do not cou- 
ple directly to external physical lines. The treatment of 
the gap equation, in a self-consistent way, is equivalent 
to a one-loop calculation. The same happens with the 
calculation of the chiral condensate and the number den- 
sity. For our purpose, we only need to obtain the Sn 
component of the propagator matrix. 

Following the construction of 40], the general Sn 
propagator is obtained in terms of the spectral density 
function (SDF): 



Sn{q) 



dk Q p(k Ql q) 
2ni k Q — q — ie 



- -n F (qo - p)p(q)- (16) 



The SDF is related to the real time propagator in Eq. 
(fl~3f through 



S (iq n ,q) = 



dk p(k ,q) 
2ni kn 



iq n 



(17) 



where the connection between both formalisms, real time 
and imaginary time, is realized through the analytic ex- 
tension of iq n — > z. What we need now is to obtain the 
SDF. In the case of free particles, the SDF can be gotten 
from the relation p — Sa(qQ + ie, q) — So(qo — ie, q)- How- 
ever, in this case where, we have a nontrivial propagator, 
we need a more general prescription to extract the SDF. 
This can be achieved by defining 



S±(q) = ± 



dz Sq(z =F ie, q) 
r ± 2iti z — qo + ie 



(18) 



where the integration path is shown in Fig. [3J With this 
definition, the SDF can be written as 



p(q) = S + (q)-S.(q). 



(19) 



In the special case of free fermions with mass M, the 
corresponding spectral density function will be 



p flee (q) = 2tt sign(q )(^ + M)8(q 2 - M 2 ). 



(20) 



By replacing pf roc in Eq. (1161) . we find the Dolan-Jackiw 
(DJ) propagator [!o[ 



Sdj(?;M) = (f + M) 



q 2 - M 2 + ie 

,2 j\/t2n 



-2wN(q )($ + M)5(q'-M 2 ) 



, (21) 




Re z 



FIG. 3. Integration path in the definition of S±. In both 
paths the integral is taken along the real axis, but the curve 
is closed through the upper half plane in F+ and through the 
lower half plane in T_. 



where the function N(qo) is defined in terms of the Fermi- 
Dirac distribution 7if(<?o) = (e q °/ T + 1) _1 and reads 

N(q ) = 9{q )n F {q - p) + 6(-q )n F (p - q ). (22) 

Note that we make a distinction here between the usual 
case where a particle is propagating freely and a more 
general situation where we have dressed propagators. 

Through this procedure, the SDF can be obtained even 
if the regulator has a cut in the real axis. Consequently, 
the Euclidean gap equation, written in Eq. (|11[) , will take 
another form in the real time formalism by introducing 
the following replacements: 



-no, 



T 



S E (q, 



(In 

E 

n 

q n ) iSu(q). 



dqo 
2tt' 



(23) 



Now, we need to construct the dressed Sn propaga- 
tors. Notice that, in the real time formalism, we are 
able to separate immediately the finite-temperature- and 
chemical-potential-dependent terms from the T, p, = 
contributions. 



IV. DRESSED PROPAGATORS FOR THE 
GAUSSIAN REGULATOR 

Following the procedure described above, we will start 
from the deconfined region, where there are two simple 
real poles. The corresponding Sn propagator is 

SH c (q) = Z{M 1 )S^{q; Mi) + Z(M 2 )S m (q; M 2 ), (24) 
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Sdj(M) being the DJ propagator with mass M and 
where Z{M) is the field-strength renormalization con- 
stant, defined as 



Z(M)= [d M i{M 2 -Y?{-M 2 )}} 
= [1 - 2M(m-M)]~ 1 . 



(25) 



In the last step, we used the pole relation M = £(— M 2 ). 
Now, this propagator describes the quasiparticles in- 
volved, and the nonlocal interaction is enclosed in the 
effective masses and the field-strength renormalization 
constants. 

The critical case that separates the real from the com- 
plex pole regions is a little bit different from the case we 
already discussed. As we can see from Fig. [51 the two real 
poles converge into a single real pole at E(0) = E c . How- 
ever, this particular case corresponds to a second-order 
pole. In addition, the new relation d q i [q 2 — S 2 (— q 2 )] = 
must also be satisfied. Following the Cauchy theorem, 
by expanding in a Laurent series and using the general 
expression that provides SDF, we find the critical prop- 
agator 



Sl?{q)=[Z c + Z' c 







S m (q;M c ), (26) 



with the constants 



Z c = 
Z' = 



2M C (4M C - 3m) 
3(2Af c - to) 2 ' 
4M 2 (M C - to) 
2M r - to ' 



(27) 
(28) 



and where M c was already defined in Eq. (|15l) . 

In both cases described above (deconfined and criti- 
cal), the finite temperature and chemical potential ef- 
fects appear separated from the T, /i = terms. Notice 
that the finite temperature and chemical potential terms 
are on mass shell due to the delta function in the prop- 
agators. This fact does not imply any difficulty when 
integrating in qo- The case of complex poles needs more 
attention. The spectral density function in this case is 
not proportional to a delta function, but it is given by a 
Breit-Wigner distribution: 



Pconf(g) 

with 



M r (A+ + A_) - i{q 2 - M 2 ){A + - A_) 



(q 2 -M 2 f+M 2 T 2 



(29) 



I ,'/,;= 2 ^!±L [q {i + M ± ) - 7 o(g 2 -A4)] ■ 



and where 



q 2 + Ml 



M 4 



A/ 2 ± iM T 



(30) 



(31) 



are the first complex solutions of the equation M — 
£(-M 2 ). Z{M) was already defined in Eq. gSJ. The 



gconf p r0 p a g a tor can be easily calculated from Eq. (fTBT) , 
but we will skip this in order to avoid long expressions. 
In fact, the full 5 xl nf will not be necessary for the rest of 
this article. 

In what follows, we are interested to separate the ther- 
mal and density effects. Then, we can write the propa- 
gator as 



Sn(q) = S (q) + S(q; T, ,u), 



(32) 



where S(q; 0, 0) = 0. In the case of the deconfined and 
critical propagators, the S term can be easily identified 
since it is proportional to N(qo). 

For the confined phase, we need to find a procedure to 
calculate integrals in the complex energy plane that in- 
volve the confined propagator, keeping in mind that when 
T and /i vanish we have to recover the Sq propagator in- 
side the integral. The finite temperature and chemical 
potential contribution must be an analytic function. We 
can see from the general expression of the Sn propagator 
that the last term in Eq. (|16l) contains all the thermal 
and density information. Nevertheless, it does not van- 
ish when T, n —¥ 0, but np(qo — fx) — > 0(—qo), which is 
not an analytic function. However, we can integrate the 
last term in Eq. (|16l) . which is analytic, proceeding then 
to remove constants independent of T and \x. Such con- 
stants then can be joined into the rest of the integral. 
This can be summarized by writing 



"Sconf (<?) = ~n F (q - fJ-)p con {(q) + C6 4 {q), 



(33) 



where C is a divergent factor that will be fixed after 
integration in order to produce a vanishing function when 
T, fj, — > 0, being absorbed into the T, \i = contribution. 
Now we proceed to calculate these integrals in order to 
obtain the gap equation, the thermodynamic potential, 
and all the other relevant quantities. 



V. GAP EQUATION 

The interesting quantities we want to calculate, such as 
the thermodynamical potential, chiral condensate, num- 
ber density, susceptibilities, etc., require momentum in- 
tegration. In particular, the thermodynamical potential 
at zero temperature and chemical potential has to be 
calculated in the Euclidean formulation. However, the 
thermal and density corrections can be handled in the 
frame of the real time formalism. 

The Gaussian regulator is constructed with the aim 
of regularizing UV divergences in the Euclidean formu- 
lation. Nevertheless, this regulator induces a divergence 
in Minkowski space when integrating in qo- The finite 
temperature and chemical potential part of the dressed 
propagators S, for the deconfined and critical case, in- 
volves a delta function which has the effect of evaluating 
the regulator at the effective masses r(— q 2 ) — > r(—M 2 ). 
The zero temperature and chemical potential term in the 
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FIG. 4. Integration path for the thermal part of the gap 
equation in the confined phase. The crosses represent the 
poles of the Fermi-Dirac distribution npfgo ± m)- 

dressed propagator So can be Wick rotated, turning back 
to Euclidean space, where the regulator produces finite 
integrals. 



Notice that the a variable and the different mass terms 
are related through M — £(— M 2 ) or, more explicitly, 

a = (M-m)e- M ^ A \ (38) 

where M stands for M^a in the deconfined phase, and 
M c in the critical phase. In the confined phase, M stands 
for all the complex poles in the confined phase, M± in Eq. 
(131 1) being the numerically relevant ones. 

For the deconfined and critical phases, the integration 
can be done in a straightforward way. However, as we 
mentioned in the last section, the situation is not as sim- 
ple when we try to integrate the thermal part of the prop- 
agator given in Eq. (|33p . In Minkowski space, the regu- 
lator in the gap equation diverges in some regions of the 
complex go plane. However, we will see that finally all 
the divergent terms do not depend on temperature and 
chemical potential and, therefore, can be removed, as will 
be explained in the last section. Replacing S from Eq. 
(1551 into Eq. J57]), we get 



We start first with the gap equation, which gives us the 
value of a. Following the previous section, we apply the 
replacements described in Eq. (p?3"|) to the gap equation 
in (fTTj) . obtaining 



a = 2N r 



**t2(-<?)trS n (q). (34) 



Now, expressing the dressed propagator in the form of 
Eq. (f52j) . the gap equation turns out to be 



MF 



da 



(35) 



where turning back to Euclidean space, the T, /i = 
contribution gives 



9o(<r) 



a N, 



G 



r 2 (p 2 )£(p 2 ) 



■2 I d PP 3 ^ ^TTZaT ' ( 36 ) 



p 2 + T l 2 {p 2 ) , 



where p = yq^ and the angular integral has already 
been done. The finite T, fi contribution in Minkowski 
space gives 



g(a,T,f,) = -2N c 



d 4 q 



r 2 (-q 2 )tYS(q;T,p). (37) 



d 4 q 



(2tt) 4 
xtrp conf (g) 



r 2 (-q 2 ) [n F (q Q - fi) + n F (q + fi)] 
2N C C 



(2nY 



(39) 



where, in the last equation, we have used the fact that 
trpconf is an odd function of qo. Now we integrate along 
the path shown in Fig. 2] If we consider only the two 
first poles of the propagator, M±, the poles inside the 

closed path will be J 'q 2 + ~M\ and -\jq 2 + M 2 . The 

same results can be obtained if we integrate in the lower 
half plane. The integrand vanishes along the upper line 
(Im go - ► oo), and the contribution from the lines sur- 
rounding the poles from the Fermi-Dirac factor (crosses) 
cancel each other. The sum of the left and right straight 
lines (Re qo — > ^foo) leaves a divergent contribution to 
the integral which, however, is independent of tempera- 
ture and chemical potential. Therefore, this term, which 
is independent of T and /i, is canceled with the constant 
factor C in Eq. ([55]). 

The deconfined and critical gap equations are obtained 
immediately, due to the presence of the delta function in 
the propagators. The resulting expressions for g are 



5(a,7» = ^]Tz(Af) f 
g CTit {a,T,fi) = IZ C + Z C — 



dkk 7 
~2E 



r 2 (-M 2 )M [n F {E - fi) + n F {E + fi)} , 



dkk 2 
~2~E r 



r 2 (-M 2 )M c [n F (E c - fi) + n F (E c + fi)] , 



(40) 
(41) 
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FIG. 5. Mi as a function of the chemical potential for dif- 
ferent values of the temperature. The end point, where the 
first-order phase transition turns to a crossover, occurs about 
{T,fx) « (70, 180)MeV. 



where k = \q\, E c = ^k 2 + M 2 7 E = + M 2 , and 
where the sum in the first equation is performed over 
all the mass poles involved in the confined phase or the 
deconfined phase. For the case of a complex mass term, 
the gap equation is real and can be easily written in terms 
of real components. 

As an example, we calculate the behavior of M\ in 
terms of the temperature and chemical potential. This 
mass term will produce the main contribution to the dy- 
namics in the deconfined phase because it is the lower 
one. The other mass term M2 turns out to be relevant 
only for values near the critical mass M c . Figure [5] shows 
the evolution of the mass Mi as a function of the chemi- 
cal potential for different values of the temperature. The 
transition observed at low temperature is obtained by an- 
alyzing the minimum of the thermodynamical potential 
as a function of a. In the next section, we will show how 
to calculate the thermodynamical potential in the real 
time formalism through the gap equation. 



VI. THERMODYNAMICAL POTENTIAL 

The inverse propagator in the real time formalism does 
not carry information about temperature and density 
[40| . This can be seen in the case of a free fermion, whose 
thermal propagator is the Dolan- Jackiw one in Eq. (PHI) , 
where the thermal information is enclosed in the on mass 
shell term. However, the inverse of the full propagator, 
~ M), does not have any information on tempera- 
ture. Nevertheless, once we have obtained the expression 
for the gap equation, we can reintegrate it, getting 



dn 



MF 



da 



da + const, 



(42) 



where the constant is independent of a and has to be cho- 
sen to regularize the thermodynamical potential. In this 
way, as we did with the gap equation, we can separate 
the zero and the nonzero contributions of the tempera- 
ture and chemical potential by defining 



fWO; T, n) = n (cr) + Q(a; T, fi), 



(43) 



where the temperature and chemical potential contribu- 
tion part obeys the relation fi(<7;0, 0) = 0. The term 
independent of the thermodynamical variables £Iq can be 
obtained directly by integrating go in Eq. (|36p. giving, 
as a result, 



n (<7) 



a 

2G 



Nc_ 
2tt 2 



dpp 3 \n[p 2 + S 2 (p 2 )]. (44) 



The finite temperature and chemical potential expres- 
sions for the gap equation, however, are not defined in 
terms of a, but in terms of the effective masses and the 
decay constant. This relation is given in Eq. (|38|) . ob- 
taining 



da = Z- 1 (M)e" M2 / A2 (iM, 



(45) 



where we have used Eq. (|25p . Proceeding with the 
integration, the temperature- and chemical-potential- 
dependent contribution to the thermodynamical poten- 
tial will be 



6(*;r,M) = -^E 



M 



dkk 2 T 



In (l + e-^-")/ T ) +/i 



(46) 



where the sum, like in Eq. (|40|) . is performed for all the 
poles involved. Both phases will be continuously con- 
nected at M = M c . This is the familiar expression for 
the thermodynamical potential. 

The thermodynamical potential with complex masses 
can be described in terms of real components. Consider- 
ing only the numerically relevant poles M±, the energy 
term can be written as 



E± = uj ±i 



MpTp 
2u ' 



with 



k 2 + M 2 + ^{k 2 + M 2 ) 2 + M 2 T 2 



(47) 



(48) 



With these definitions, we can write the finite T, \i contri- 
bution to the thermodynamical potential in the confined 
phase as 



6 (~ t \ ANc 

S iconic; 1 , fl) = 5- 



dkk 2 T 



x [ln(l + B_)+ln(l + £+)], (49) 
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with the functions B± defined as 

/M r 



2 cos I ) e ^ ± ^' T + e - 2 ^)/ T 

V 2Tw 



(50) 



Our construction was compared with |23| . obtaining the 
same results. The cosine term in Eq. (150]) can produce 
an unstable thermodynamical potential [111 ] . However, 
the inclusion of the Polyakov loop should fix this problem 
[25l l3lj . Note that now, we can also obtain other relevant 
quantities directly from the thermodynamical potential 
already calculated in the real time formalism. The quark 
number density function is obtained by taking the deriva- 
tive with respect to the chemical potential: n = — 9il/<9/i. 
The chiral condensate is calculated through the deriva- 
tive of the regulated thermodynamic potential with re- 
spect to the mass: (qq) = dfl r /dm. Here, the current 
mass is related to the effective mass terms through Eq. 
(|38|) . and the partial derivative is done by taking a and 
A as constants. 



tion [33, |40J , we generalize the procedure to get the spec- 
tral density function. This generalization allows us to 
deal with dressed propagators whose analytical struc- 
ture includes complex poles. With this, we obtain differ- 
ent propagators if the quasiparticles are deconfined (real 
poles), confined (complex poles), and also the critical case 
which separates both regimes. Once we have obtained the 
real time propagators, we find the gap equation which 
provides the value of a and, consequently, the value of 
the quasiparticles effective mass terms. The thermody- 
namical potential is obtained from the gap equation, and 
it is given in terms of simple expressions as a function 
of the effective mass terms. We verified that our results 
coincide with the ones reported in [23| . 

This procedure gives an intuitive phenomenological de- 
scription in terms of quasiparticles. Its generalization to 
other regulators is, in principle, straightforward. In fact, 
our prescription for obtaining the spectral density func- 
tion is well-defined, even if the propagator presents a cut 
along the real axis, which is the case for some type of 
Lorentzian regulator with fractional power. 



VII. CONCLUSIONS 
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